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O ■ Abstract 

In this article, we set up a functional setting for mean-field electronic struc- 

Cd i ture models of Hartree-Fock or Kohn-Sham types for disordered crystals. The 

electrons are quantum particles and the nuclei are classical point-like particles 

whose positions and charges are random. We prove the existence of a minimizer 

of the energy per unit volume and the uniqueness of the ground state density 

of such disordered crystals, for the reduced Hartree-Fock model (rHF). We con- 

T^ ' sider both (short-range) Yukawa and (long-range) Coulomb interactions. In the 

^pH. former case, we prove in addition that the rHF ground state density matrix 

r^ ' satisfies a self-consistent equation, and that our model for disordered crystals is 

the thermodynamic limit of the supercell model. 
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1. Introduction 

The modeling and simulation of the electronic structure of crystals is one of 
the main challenges in solid state physics and materials science. Indeed, a crystal 
contains an extremely large number (in fact an infinite number in mathematical 
models) of quantum particles interacting through long-range Coulomb forces. 
This complicates dramatically the mathematical analysis of such systems. 

Finite size molecular systems containing no heavy atoms can be accurately 
described by the A''-body Schrodinger equation, or its relativistic corrections. 
Because of its very high complexity, this equation is often approximated by 
nonlinear models which are more amenable to numerical simulations. On the 
other hand, no such reference model is available for infinite molecular systems 
such as crystals. For this reason, in solid state physics and material sciences, 
the electronic structure of crystals is often described by linear empirical models 
on the one hand, and mean-field models of Hartree-Fock or Kohn-Sham types 
on the other hand. 

In linear empirical models, the electrons in the crystal are seen as non- 
interacting particles in an effective potential VeS, so that their behavior is com- 
pletely characterized by the effective Hamiltonian 

a self-adjoint operator on L^(K.''). Here d is the space dimension which is (i = 3 
for usual crystals. The cases d ~ 1 and d = 2 are also of interest since linear 
polymers and crystalline surfaces behave, in some respects, as one- and two- 
dimensional systems, respectively. Throughout this article, we adopt the system 
of atomic units in which h ~ 1, m^ = 1, e = 1 and 47reo = 1, where h is the 
reduced Planck constant, me the mass of the electron, e the elementary charge, 
and £q the dielectric permittivity of the vacuum. For the sake of simplicity, 
we work with spinless electrons, but our arguments can be straightforwardly 
extended to models with spin. 

When the system under study is a perfect crystal, the effective potential 
Vcs is an 7?,-periodic function Vpor, where 7?, is a discrete lattice of R'^, and 
the effective Hamiltonian is then a periodic Schrodinger operator on L^(E'*), 
H = gpcr = — ^A -|- Vper- The spectral properties of such operators are well- 
known |3(|]. Under some appropriate integrability conditions on Vpcr, it follows 
from Bloch theory that the spectrum of -ffpor is purely absolutely continuous 
and composed of a countable number of (possibly overlapping) bands. 



It is possible to describe local defects in such effective linear models. Displac- 
ing or changing the charge of a finite number of nuclei corresponds to adding a 
potential W to Vpor- Because such perturbations are local, the potential W de- 
cays at infinity and therefore the effective Hamiltonian fldefect = — ^A-l-ypor+W^ 
has the same essential spectrum as the unperturbed Hamiltonian Hpa. On the 
other hand, ifjefect niay possess discrete eigenvalues below its essential spec- 
trum, or lying in spectral gaps. They correspond to bound states of electrons 
in the presence of the local defects. 

Doped semiconductors and alloys are examples of disordered crystals, which 
are perturbed in a non-local fashion. Such systems can be adequately modeled 
by random Schrodinger operators |7|, |35| . One famous example is the continuous 
Anderson model 

H^ = --A + V^ with V^{x) = J2li'^i^)xix-k), 

ken 

where, typically, x ^ C^IR"^) and the g/^'s are i.i.d. random variables. Here, 
only the charges are changed but it is possible to also account for stochastic 
displacements. The study of the spectral properties of ergodic Schrodinger 
operators is a very active research topic (see e.g. [16| and the references therein). 

In linear empirical models, the interactions between electrons are neglected 
(apart from the implicit interaction originating from the Pauli principle pre- 
venting two electrons from being in the same quantum state). Taking these 
interactions into account is however a necessity for a proper physical descrip- 
tion of these systems. One main difficulty is then that the Coulomb interaction 
is long-range and screening becomes extremely important to explain the macro- 
scopic stability of such systems. Understanding screening effects in a precise 
manner is a difficult mathematical question. 

As already mentioned above, there is no well-defined many-body Schrodinger 
equation for crystals. The only available way to rigorously derive models for 
interacting electrons in crystals is to use a thermodynamic limit procedure. The 
idea is to confine the system to a box, with suitable boundary conditions, and 
to study the limit when the size of the box grows to infinity. For stochastic 
many-body systems based on Schrodinger's equation, it is sometimes possible 
to show that the limit exists. In [37|, Veniaminov has first considered a many- 
body quantum system with short range interactions. Short after, the existence 
of the limit for a crystal made of quantum electrons and stochastic nuclei in- 
teracting through Coulomb forces was shown in [3|, by Blanc and the third 
author of this article. In these two works dealing with the true many-body 
Schrodinger equation, the value of the thermodynamic limit is not known. For 
Thomas-Fermi-type models, Blanc, Le Bris and Lions were able to identify the 
thermodynamic limit and to study its properties [2|. Unfortunately, Thomas- 
Fermi theory is not able to reproduce important physical properties of stochastic 
quantum crystals, like the Anderson localization under weak disorder. 

The purpose of the present work is to propose and study a mean-field 
(Hartree-Fock type) model which can be obtained from a thermodynamic limit 
procedure, for an infinite, randomly perturbed, interacting quantum crystal. 
This model is not as precise as the many-body Schrodinger equation, but it is 
still much richer than Thomas-Fermi type theories. In particular, it seems ade- 
quate for the description of Anderson localization in infinite interacting systems. 



More specifically, we consider a random nuclear charge ^{uj,x) > 0. For 
simplicity we do not consider point-like charges, and we assume that //(a;, •) e 
Ll^^{M.'^) almost surely. Also we are interested in describing random perturba- 
tions which have some space invariance, and we make the assumption that they 
are the same when the system is translated by any vector of the underlying peri- 
odic lattice TZ. We assume that the group TZ acts on the probability space in an 
ergodic fashion and we always make the assumption that fi is stationary^ which 
means fi{Tk{uj)^ x) = fi{uj, x + fc), where t = {Tk)k<£'R is the ergodic group action 
on the probability space. A typical example is given by a lattice TZ with one 
nucleus per unit cell, whose charge and position are perturbed by i.i.d. random 
variables, 

M(w,a;) = ^ qk{i^)x{x- k-ijk{^)). 
ken 

Similarly, the state of the electrons in the crystal is modelled by a one-particle 
density matrix |23| . that is, a random family of operators 7(0;) : Li^iM?) — >■ 
L/^^E?) such that < 7(0;) < 1 almost surely. It is also assumed that 7 is 
stationary in the sense that its kernel satisfies ^{Tk{uj), x, y) = 7(0;, x -f- fc, y + fc) 
for all k £ TZ. These concepts will be recalled later in Section [Ol 

For any such electronic state 7 we define in Section U] the corresponding 
reduced-Hartree-Fock (rHF) energy, in the field induced by the nuclear charge 
fj,. This energy is just the sum of the kinetic energy per unit volume of 7 and the 
potential energy per unit volume of 7 and fi. The rHF model is obtained from 
the generalized Hartree-Fock model [2J, ll| by removing the exchange term [3J| . 
Alternatively, it can be seen as an extended Kohn-Sham model [l^l with no 
exchange-correlation . 

Defining the rHF energy properly requires to introduce several tools, which 
is the purpose of Sections [5] and |21 We start by defining the average number 
of particles and the kinetic energy per unit volume for ergodic density matrices 
and we show useful inequalities. In particular we derive Hoffmann-Ostenhof |l7| 



and Lieb-Thirring inequalities |26l . l27j for ergodic density matrices, which are 
very important estimates that we use all the time. Loosely speaking, they can 
respectively be stated as follows: 



Average kinetic energy per unit vol. of 7 > El / |Vy^| 



and 



n5/3 



Average kinetic energy per unit vol. of 7 > KE I / pi? 

where Q is the unit cell, p^ is the density of the state 7 and X is a constant 
independent of 7. 

In Section [21 we discuss Poisson's equation 

- Ay = Anp (1) 

for stationary functions p{ijj, x), where A is the Laplace operator with respect to 
the x- variable, and we explain that the situation is much more complicated than 
in the periodic case. In particular, the neutrality condition E( J„ p) = on the 
charge density appearing on the right side of ([Ij is necessary but in general not 
sufficient to find a stationary solution V. When E{J„ p^) < cx3 and IE(/q p) = 0, 



it is possible to give a necessary and sufficient condition for the existence of 
a stationary solution T^ to ([1]) such that IE(/n V^) < oo. In words, p should 
be in the range of the "stationary Laplacian" which is a particular self-adjoint 
extension of —A on L^{i^ x Q) with "stationary boundary conditions". 

Understanding Poisson's equation ((ij for general stochastic charge densities 
p is an important and interesting problem in itself. In order to define the 
associated Coulomb energy per unit volume, we adopt here a simple strategy 
and take the limit ?7i — >■ of the Yukawa energy. This means we consider the 
regularized equation 

^AVrn+m^Vm=47Tp (2) 

and we define the Coulomb energy as the limit of E ( /„ Vmp) when m ^- 0. We 

then give in Section [3] several properties of this energy. 

After these preliminaries, we are able to properly define and study the re- 
duced Hartree-Fock energy for stochastic crystals in Sectional In particular we 
prove the existence of a minimizer 7 of this energy and the uniqueness of the 
ground state density Pj. In the Yukawa case to > 0, we also show that the 
minimizers solve a self-consistent equation of the form 

-AV„-, + m^Vra = 47r(p-, - p) . 

The mean-field operator 

H„i = —-^^A + V,n 

is a random Schrodinger operator describing the collective behavior of the elec- 
trons in the system. Studying its spectral properties would allow to understand 
localization properties in the interacting stochastic crystal. 

In Section [SJ we finally prove that, in the Yukawa case, our model is actu- 
ally the thermodynamic limit of the supercell reduced Hartree-Fock theory (the 
system is confined to a box with periodic boundary conditions). This justifies 
our theory with Yukawa interactions. For Coulomb forces, our proof does not 
apply because of some missing screening estimates. We make more comments 
about this later in Section [51 

Let us end this introduction by mentioning that our theory is rather general 
and it actually works for any reasonable interaction potential which decays fast 
enough at infinity. We concentrate on the Yukawa interaction because of the 
limit TO, — > which corresponds to the more physical Coulomb case and which 
we study as well in this paper. Note that we consider here the action of a 
discrete group on Q because we have in mind the case of a randomly perturbed 
crystal. Our approach can also be applied to the case when the group acting on 



fl is R (amorphous material). We refer to [18| for details 
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2. Electronic states in disordered crystals 

In mean-field models (such as Hartree-Fock or Kohn-Sham), the state of the 
electrons is described by a self-adjoint operator 7 acting on L^(R'^), satisfying 
< 7 < 1 in the sense of quadratic forms, and such that Tr (7) is the total 
number of electrons in the system [23|. In (infinite) crystals, we always have 
Tr (7) = -|-cx). Such an operator 7 is called a (one-particle) density matrix. The 
purpose of this section is to recall the main properties of electronic states in 
a class of random media, satisfying an appropriate invariance property called 
stationarity. 

2.1. Basic definitions and properties 

Throughout this paper, d will denote the space dimension. We will later 
focus on the cases where d £ {1, 2, 3}, but we keep d arbitrary in this section. 
We restrict ourselves to the cubic lattice group 7?, = Z'^ to simplify the notations; 
general discrete subgroups TZ can be tackled similarly without any additional 
difficulty. We consider a probability space (17, J^, P) and an ergodic group action 
T oilA onVl. We recall that r is called ergodic if it is measure preserving and if 
for any A e J" satisfying Tk{A) = A for aU k e Z'\ it holds that P(A) e {0, 1}. 

Example 2.1 (i.i.d. charges). A typical probability space we have in mind is 
the one arising from a random distribution of particles of charges qi and q2 on 
the sites of the lattice H^ with probabilities pi and 1— pi. The probability space 
is then given by il = {qi, (72}^ (md P = p®^ where p = PiSq^ + (1 — Pi)Sq2 . In 
this case, the group action is Tk{uj) = w.+fc. 

The ergodic theorem |36l . Theorem 6.1, Theorem 6.4], which will be exten- 
sively used in the sequel, can be stated as follows: 

Theorem 2.2 (Ergodic theorem). If t is an ergodic group action of U^ on il 
and X G U'{Vt), with 1 <p < 00, then, 

hm — J— ^ ^ X{Tk{Lo))=nX). 

"^°° 2n+l ,V 1 

almost surely and in L^ (f2) . 

A measurable function / : il x M** — > C is called stationary if 

Vfc G Z'', f{Tk{uj),x) = f{u},k + x), a.s. and a.e. 

We will make use of the families of stationary function spaces 

LP (L«) = {/ e LP {n, Ll^ {R")) I / is stationary} , 

and 

HT = {f (^L^ {n, Hl2c (K"^)) I / is stationary} , 

and resort, for convenience, to the shorthand notation L^. = L^ [L^). Endowed 
with the norms 



E 



LUL") =^niillL<!(Q) 



(ll/lli. 



and the scalar products 

(/,g)i. =E((/,5)i2(Q)), (/,5>Hr =E((/.3>H'"(Q)) 

where 

^ [ 2' 2 

denotes the semi-open unit cube, the spaces if (i*) are Banach spaces and the 
spaces Lg and H^ are Hilbert spaces. 

We denote by "H = i^ (R'') the space of complex valued, square integrable 
functions, equipped with its usual scalar product (•, ■)l'^. We also denote by 

• B the space of the bounded linear operators on "H, endowed with the 
operator norm |j-|j; 

• S the space of the bounded self-adjoint operators on "H; 

• &p the p'^ Schatten class on H. Recall that ©i is the space of the trace 
class operators on T-L and ©2 the space of the Hilbert-Schmidt operators 
on H. 

Let 2? be a dense linear subspace oiJi. A random operator with domain 2? is a 
map A from J7 into the set of the linear operators on H such that V C D{A{lj)) 
a.s. and such that the map uj 1— > {A{uj)x^ y) is measurable for all x ^ T) and 

Of importance to us will be the uniformly hounded random operators A which 
are such that supess^^gn ||^(<^)|| < 00. The Banach space of such operators is 
denoted by L°°{Q,,B). This is a M^*-algebra which is known to be the dual of 
L^{Vt,&i) (see, e.g., [3l|, Corollary 3.2.2]). We will often use the corresponding 
weak-* topology on L°°{Q,, B) for which An ^* A means 

E(Tr(A„B)) ^E(Tr(ylS)) 

for all B G L^(il,Si). Since L^{Q,,&i) is separable, any bounded sequence 
(An) in L°°(0,S) has a subsequence [An^) which converges weakly-* to some 
A S L°°{Q,,B). Similarly, we know that the dual of L'P{Q,,&q) is nothing else 
but LP' {Vl,&q') when 1 = l/p+l/p' = l/q+l/q' and 1 < p, g < 00. 

Let (Uk)k£i.^ be the group of unitary operators on "H defined by 

Ukfix) = fix + k), a.e., V/ en,yke Z'^. 

A random operator A (not necessarily uniformly bounded) is called ergodic or 
stationary if for any fc G Z'', I? C Uki'D) and the following equality holds 

A{TkiLj)) = UkA{io)U*k, a.s. 

One of the fundamental theorems for ergodic operators [35|, Theorem 1.2.5 p. 13] 
states that for any self-adjoint ergodic operator A, there exists a closed set 
S C R and a set f7i e J" with P(r2i) = 1, such that cr(A(a;)) == E, for aU 
Li! G ill. The set E is called the almost sure spectrum of A. 

We finally denote by S_ the space of the ergodic operators on H that are 
almost surely bounded and self-adjoint. 



2.2. Ergodic locally trace class operators 

In this section, we recall the definitions of the trace per unit volume, the 
density and the kernel of an ergodic locally trace class operator (see e.g. |4l.lll|). 
For 1 < p < oo, we denote by L^ (W^) the space of the compactly supported L^ 
functions on M''. 

Definition 2.3 (Locally trace-class operators). A random operator A is called 
locally trace class if x^X € L^{^,&i) for all x e L'^(R'^), that is, 



VxeLr(R'), IE(Tr(|xA(.)x|)) < 



We now focus on the particular case of ergodic operators, and denote by &i 
the space of the ergodic, locally trace class operators. The following character- 
ization of the positive operators of &-^ will be useful. 

Proposition 2.4 (Characterization of ergodic locally trace-class operators). 
Let A be a positive, almost surely bounded, ergodic operator. Then A is locally 
trace class if and only z/E(Tr (Iq^(-)Iq)) < oo. 

The trace per unit volume of an operator A G 6-^ is defined as 



Ti(A)=E(Tr(lQ^(.)lQ)). (4) 



The following summarizes the main properties of locally trace-class ergodic 
operators. 

Proposition 2.5 (Kernel and density). Let A £ ©j^. Then, there exists a 
unique function A{-,-,-) E L^{il, Lf^^{M.'^ x M'')), called the kernel of A, and a 
unique function pA G Ll , called the density of A, such that 

Wip e LI{R'^), {A{uj)ip) (x) ^ / A{uj,x,y)ip{y)dy a.s. and a.e. 

and 

Vxeir(R'), Tr(xA(c.)x)= / x'{x)pAiLO,x)dx a.s. (5) 

The kernel A{-, •, •) is stationary in the following sense 

A{Tk{u!),x,y) = A{uj,x + k,y + k), yk E Ij a.e. and a.s. 
Moreover, if A>0, then pA > 0. 

Note that it follows from (g]) and © that 



Tr (A) = E (^1 



PA 

The proofs of Propositions 12.41 and 12.51 are elementary; they can be read in |18j . 

The following cyclicity property is proved in jll| , based on arguments in [lO| 
(see also [18| for a detailed proof): if B is an ergodic operator in L°° {^,B) and 
A a positive operator in ©^ nL°°(fi,S), then BA and AB are in ©^ nL°°(fi,S) 
and 

Ir {BA) == Tr {AB) . (6) 



2.3. Ergodic operators with locally finite kinetic energy 

Ergodic density matrices for fermions are operators 7 S Sj^ n S_ such that 
< 7 < 1 a.s. By BirkhofF's theorem, the trace per unit volume can be 
interpreted from a physical viewpoint as the average number of particles per 
unit volume. In this section, we define and study in a similar fashion the average 
kinetic energy per unit volume. 

2.3.1. Definition 

For 1 < j < d, as usual, we denote by Pj — —idx the momentum operator 
in the j"^ direction, which is self-adjoint with D{Pj) — {^p & % \ dx ^p € %}. 
As Pj commutes with the translations, we see that for all A e &_^, the operator 
PjAPj is ergodic. The operator PjAPj is well defined and bounded on D{Pj), 
with values in D{Pj)' , where D{Pj)' is the topological dual space oi D{Pj). We 
say that the kinetic energy of A is locally finite if PjAPj G &i, and we then call 



Tri-AA):=Y,Tr{P,AP,] 



the average kinetic energy per unit volume of A. We denote by &i ^ the subspace 
of &i composed of the ergodic locally trace class operators with locally finite 
kinetic energy. 

2.3.2. Hoffmann-Ostenhof and Lieb-Thirring inequalities for ergodic operators 
For finite systems h € 61 H 5, < 7 < 1 and Tr (-A7) < 00), the 
Hoffmann-Ostenhof (iTI . |23| and Lieb-Thirring [26J, |27|, |23| inequalities provide 
useful properties of the map 7 1— > p^. In this section, we state and prove an 
equivalent of these inequalities for ergodic density matrices with locally finite 
kinetic energy. 

Proposition 2.6 (Hoffmann-Ostenhof inequality for ergodic operators). Let A 
be a positive operator in 61 1 H 5. Then 

y^ e Hi and E (J \\7^^A < Tr (-AA). 

Proof. It follows from Proposition l2.5l that ^fpA G L^. Let i? be a compact set of 
R'^ and 7? G C;?°(M'^) such that 7? = 1 on B and < 77 < 1. The operator T]A{Lu)'q 
has finite kinetic energy a.s. Therefore, the Hoffmann-Ostenhof inequality gives 



|V^PA(w)| == |V^p,,A(cu)r,| < /X! A„(a;)|V<^„(w)|2 a.s. and a.e. on B, 



where {ipn{uj))^^j^ is an orthonormal basis of eigenvectors of the compact self- 
adjoint operator r]A{u!)r] and (A„(w))„gj^ the associated eigenvalues. As 

« d 

/ ^A„H|V^„H|2=^Tr (lBP,r/A(tj)7?P,lB) a.s. 

•'^ TiSN 3 = 1 



and as for all 1 < j < d^ 1b [Pj, v] = —ilBdx^ = 0, we deduce that 

^. d d 



Therefore 

E 



Jb J _i 



As ^ has locally finite kinetic energy, we conclude that ^/pa G Hi. For B = Q, 
we obtain the stated inequality. D 



The following corollary is an obvious consequence of Proposition 12.61 and of 
the Sobolev embeddings. 

Corollary 2.7. Let A be a positive operator in &i i O S_. Then, pA G Ll [V), 
for p = +00 if d=l, pe [l,+oo) if d = 2 and 1 <p< -^ if d>3. 

The following is now the ergodic equivalent of the Lieb-Thirring inequal- 
ityS [23,123. 



Proposition 2.8 (Lieb-Thirring inequality for ergodic operators). There exists 
a constant K{d) > 0, depending only on the space dimension d > 1, such that 
for all 7 e 6-^ 1 n 5 with < 7(w) < 1 a.s., 

p^ e L7^ and K{d) E ( f p/^ j < Ir (-A7). (7) 

Proof. To prove ([7|) , we apply the Lieb-Thirring inequality in a box of side- length 
L, and then let L go to infinity. The constant K{d) can be chosen equal to the 
optimal Lieb-Thirring constant in the whole space. Let F^ ~ [— L/2,i/2) 
and let (xl)^^^. be a sequence of localizing functions in C^ (R'^), such that 
< XL < 1, XL = 1 on Fz,_i, XL = outside of Tl, and |VxL(a;)| < C. We 
first apply the Lieb-Thirring inequality to Xl7('^)xl and obtain 

K{d) pj{uj,x)'^ dx<TT {-Axli{uj)xl) a.s. 

Next, using the stationarity of p..y and the equality [Pj,xl] — —id-xXL, we get 
for any £ > 

i^(d) E f / pT") < -il±ii^ ^ E (Tr (xlP,7^,Xl)) 

+ 7^EE(Tr((^^.XL)7(5.,XL))). (8) 

For each 1 < j < rf, we have 

Tr (xlP,7HP:,Xl) < E Tr (lQP,-7(r, H)P,1q) a.s. 
feertnz'* 



10 



It follows from the ergodicity of 7 (hence of Pj-jPj) that 

d 

J2 E (Tr ixLPnP.XL)) < L^Tr {~A^) . (9) 



Besides 
Tr 



{{dx,XLh{dx,XL)) = P',{i^r){dx,XL)^ <C Pji^r), 

where we have used that Vxl is uniformly bounded. Using again the station- 
arity of p^ , we obtain 

d 

^E (Tr (dx^XL) 7 {d.,XL)) < CL'^-'Tr (7) • (10) 

Combining ([5]), ^ and pTI)l . letting L go to infinity then letting e go to 0, we 
end up with the claimed inequality. D 

2.3.3. A compactness result 

In this section we investigate the weak compactness properties of the set of 
fermionic density matrices with finite average number of particles and kinetic 
energy per unit volume 

/C := {7e 6i_in5, <7 < 1 a.s.}. (11) 

This set is a weakly-* closed convex subset of L°°{il,B). The following result 
will be very useful. 

Proposition 2.9 (Weak compactness of ergodic density matrices). Let (7„) be 
any sequence in /C. Then there exists 7 G /C and a subsequence (n^) such that 

1. In, ^* jtnL^in,B), 

2. lim^(7„J=^(7), 

d + 2 

3. Pj^^ — ^ p-y weakly inLg'' , 

4. Tr (-A7) < lim infix (-A7„). 

n— foo 

Note that, in average, there is never any loss of particles when passing to 
weak limits: Ir ("/„) tends to Ir (7) as n — > 00. On the other hand, even if 
we have p^^ -^ p^ weakly and K{J„ p-y^J — > E(J„p^), in general we do not 
have almost sure convergence and we do not expect strong convergence in Lf 
for 1 < p < 1 + 2/d. 

Example 2.10 (Weak versus strong convergence for p^^). Consider a smooth 
function if with compact support in the ball -6(0, 1/2) such that \\(p\\j^2 = 1, and 
the operator 

'- '- \ipi- + k)){ip{- + k)\. 

keZ'' 

11 



where (ajfc)^,,^^ are i.i.d. variables, uniformly distributed on [0,1]. Then we 
have 7„ G ©^ i, < 7„ < 1, 



1^ |^(. + fc))(^(. + fc)| ^nL^in,B) 



2 

fcGZ'' 



P7« 



El + sin(27ma;fe) , , ,,,2 1 v^ 1 , ,,1 



weakly in Lf /or 1 < p < 00 . We also have 



E ( / p^ I = E ( / p^ I , Vn e N 

However, since sm{2Tmujk) -^ weakly but not strongly in £''([0, 1]), we do not 
have any strong convergence for p^^ . 



Proof of Provosition \2.y[ Consider a sequence (7^) as in the statement. Since 
(7„) is bounded in L°°{il,B), there exists 7 € L°°{^,B) such that 7„ converges 
to 7 weakly-* in i°°(J7, B), up to extraction of a subsequence (denoted the same 
for simphcity). Recall that 7„ ^* 7 means 

lim E (Tr (A-fn)) = E (Tr (A7)) 

for all j4 e L^{il,&i). Using for instance A — Y\f){g\ for some fixed f,g £ 
% = L^(]R'') and some fixed Y e L^{il), we find in particular that 

yf,gen,yYeL\n), E(y(g,7/))= lim E(r(5,7„/)). (12) 

Hence, {gj-fnf) converges to {g,jf) weakly—* in L°°{fl). Using this, it is easy 
to verify that 7 is ergodic and satisfies 7* = 7, < 7 < 1 a.s. 

Let now {fk)k>i be any orthonormal basis of LF'{Q) where we recall that 
Q is the unit cell. Using that E((/fc,7„/fc)) — > E ((//;, 7/^)) for each fc > 1 as 
71 — > 00, and Fatou's lemma in ^'"'^(N), we obtain 



E(Tr(lQ7lQ)) = VE((/fe,7/,)) < liminf E V (A,7„/,) 

^ — ^ ri— >-oo \ ■'- — ^ 

fc>l yfe>l 

= liminfE(Tr(]lQ7„lQ)). 

n— ^00 

By Proposition 12.41 we conclude that 7 e @_i . The same argument can be 
employed to show that 7 G &_i j^, assuming this time that each /^ is in Hq{Q). 
Then we have for each k 

lim E((/fe,P,7„P,/fe))= lim E (((P,-/fe), 7„(P,/fe))) = E (((P,A), 7(P,/fe))) , 
by (fT^ and with Pj — —idx- By Fatou's Lemma in £^{N) we see that 

d 

Tt (-A7) = V VE(((P_,-A.),7(P,-/,))) < liminf^ (-A7n) ■ 

^ — ^ ^ — ^ n— >-oo 
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Let us now prove that Tr (7,1) indeed converges to It (7). We consider a 
smooth function x in C'^{W^). The sequence (7„) being bounded in &_^ ^, there 
exists a constant C £ M+ such that for ah 71 G N and 1 < j < d, 

E (Tr (x7nX)) + E (Tr {xFjlnP^x)) + E (Tr (9,^x) 7n (a.^x)) < C. 
Using again the relation [-P,,x] = —idxX^ '^^ obtain 

E(Tr(P,x7nXP,))<4C, 
hence 

E (Tr ((1 - A)^ xinX (1 - A)^)) = E (Tr {xinX)) + E (Tr (-A(x7nX))) 

d 
= E (Tr (x7nX)) + 51 E (Tr (P,x7nX-P,)) 

< (1 + Ad)C. 

This proves that (1 — A) ' X7nX (1 ~ A) ' is bounded in L^(p,, &i) or, equiv- 
alently, that (1 — A) ' Xy/ln is bounded in L^{^, ©2)- From this we infer that 

(1 - A)^ x7nX = {(1 - A)^ x7nX (1 - A)^} (1 - A)"^ 

is bounded in i^(57, 61), since (1 — A) ' is a bounded operator. Similarly, 
we can write 

(1 - A)^ X7„X = {(1 - A)^ xVt;^} Vi7.X, 

which is now bounded in Lp^{^^ ©2), since k/TnX ^ ^ due to the assumption 
that < 7„ < 1. We conclude that (1 — A) ' X7nX is bounded in i^(0, ©1) n 
L?{yt, ©2), hence in LP{VL, ©p) for all 1 < p < 2, by interpolation. In particular, 

(1 - A)^ X7«X ^ (1 - A)^ X7X weakly in LP(f7, ©p) for aU 1 < p < 2. (13) 

That the limit can only be (1 — A) ' X7X follows for instance from P^ with 
functions /, .g e H'^{W^). 

We consider now a fixed function Y e L°°{^) and write 

E (FTr (x7nX)) - E (rTr f (1 - A)^/' X7«x1b (1 - A)"^/" 



where i? is a large enough ball containing the support of x- By the Kato-Seiler- 
Simon inequality |33l . Thm. 4.1], 

Vp > 2, ||/(x•)5(-^V)||e^ < {27:)-^/^ ||/ILp(e-) Ml^^^) , (14) 

we have 



Ib(I-A)"^/^ < {2t:)-'^Ip\B\^ ( f 



dp 

l + d 



1 + d 



1+rf WR'' (1 + IpP) 
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hence Is (1 - A)"^^^ e Gi+d- Thus Fls (1 - A)"^/^ e L°°(f],6i+d) c 
L^+''(J7, &i+d)- Since 1 < 1 + l/d < 2 we obtain by the weak convergence P^ 

in Li+Vrf(f], 6,^,/,) = L^+d^n, ei+d)', 

J[|m E(rTr(x7nX))=E(rTr ((1 - A)^/' x7X1b (1 - A)-^/')) =E(rTr(x7X)) 
We can reformulate this into 



n-^oo 



hm E r / x>7„ ]=^[Y X>7 . (15) 



for all Y G i°°(r2) and aU x e C;?°(R''). 

As Tr (— A7„) is bounded, we infer from the Lieb-Thirring inequality for 
ergodic operators (Proposition 12. 8p that (p^^ ) is bounded in L^ . We can 
therefore extract a subsequence which weakly converges in Ls to some p G 

Ls . Since the space spanned by the functions of the form Yx^ with Y e 
L°°{Q) and X e C^iQ) is dense in L^+'^/^{n x Q), we deduce from dH]) that 
p^ = p. Now. using that Ig S L^+'^/^(r2 x Q), we finally obtain the claimed 
convergence 

hm S (7„) = lim E ( / p^„ ) = E f / p^ ) = TV (7) . (16) 

This concludes the proof of the proposition. D 

2.3.4- Spectral projections of ergodic Schrodinger operators 

The following result provides a control of the average number of particles 
and kinetic energy per unit volume of the spectral projections of an ergodic 
Schrodinger operator, in terms of the negative component V- = max(— V,0) 
of the external potential. We will use it later in Section 14.41 to prove that the 
ground state density matrix of the reduced Hartree-Fock model with Yukawa 
potential is solution to a self-consistent equation. 

Proposition 2.11 (Spectral projections are in &_^ ^). Let V € L^ he such that 
the operator H = —A + V is essentially self-adjoint on C^{M.'^) and V- G 
Ls . Denote by Px ~ l(-cc>.A) (H) the spectral projection of H corresponding 
to filling all the energy levels below X. Then, Px G &i i for any A G K and there 
is a constant C > (depending only on d > 1) such that 

d 

Tr{Px)<C (e( f {V-X)^)Y^\ (17) 



Tr{~APx)<CE(f{V-X)Jj. (18) 

The estimate (flT)) on Tr (Px) is probably not optimal but it is sufficient for 
our purposes. 
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Proof. Let us first prove that P\ <E &i i under the assumption that VL G 
L°° (p, X M''). By the Feynman-Kac formula [32, Theorem 6.2 p. 51], we have 
for alH > 

p^-n-A+v-, < — ,,„ a.s. (19) 

(47rt)^ 

Using then the inequality 1(_oo,a] (2;) < e^*'^^'^^ for all A G M and all t > 0, 
as well as the fact that V is uniformly bounded from below, we deduce that 
P\ G &i and pp^ G L°° (jl x R''). Likewise, using the inequality x1(_oo,a] {x) < 
AlA>oe"*^^"^\ we obtain that HPx G ©i, hence that Px G ©ii- 

Now that we know that Pa G &i j^, we can derive bounds which only de- 
pend on ||(y — A)_|Li+d/2. The general case will then follow from a simple 
approximation argument. We start by noting that 

0<Tr ((-A + V^-A)_) = -Ii ((-A + 1/ - A) Pa) 

< ~C \\ppA%2 + WppJ ^\\{v~x)4 ^ 

T — 7r~ r d '' '' T 2 

< C||(l^-A)„||'tL> 

where we have used the Lieb-Thirring inequality ([7]) for ergodic operators. 
Therefore 

Il{Px) = \\ppJl^<C\\ppJ ^ <c||(y-A)_||^+, . (20) 

As < Tr ((-A + V - A)_), we obtain 

S(-APA)<C||(y-A)_||T+, . (21) 

This concludes the proof in the case of bounded below potentials. In the general 
case we consider the sequences of cutoff potentials Vn ~ max {V, ^n} and corre- 
sponding operators Hn = — A -I- y„ and show that for any bounded continuous 
function /, the operator / (i?„) converges to / (H) in the strong operator topol- 
ogy a.s. We conclude the proof using an appropriate approximation of 1(_oo.a] 



by bounded continuous functions (see [18| for details). D 



We can now use the previous theorem to deduce a useful variational char- 
acterization of the spectral projection P\ among all ergodic fermionic density 
matrices 7 G /C having a locally finite kinetic energy. 

Proposition 2.12 (Variational characterization of spectral projections). As- 
sume that V is as in Provosition \2.1I\ and denote again P\ := 1(_oo.a)(^) with 
H — —A + V . For every A G R, the minimization problem 

inijTr{-A^)+E(Jvp^)^XTT{j)\ (22) 

admits as unique minimizers the operators of the form J — P\ + 6 where < 
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Note that K{J„V p^) is well defined in (—00, +00] since V- £ Ls by 

assumption, whereas p^ E Lg by the Lieb-Thirring inequality ([7]). 

Proof. When 7 is smooth enough (— A7 G ©j for example) and V £ Lf, we 
can write 

Ty (-A(7 - Pa)) +¥.(j V {p^- pp,)\ - A^ (7 - Px) 

- lL ((-A + V- A)(7 - Pa)) > :n (I - A + y - A|(7 - Pa)^) . 

In the last estimate we have used the cyclicity property © and the fact that 

Pa^(7 - Pa)Pa^ -Px{i- Px)Px > (7 - Px)\ 

which turns out to be equivalent to < 7 < 1. A simple approximation argu- 
ment now shows that the inequality 

TY{~A-f)+E(f Fp^Va^(7)>Ti(-APa)+e('/' Vpv^-XTl(Px) 

+ Tt (\H - X\'/'{j - P^f\H ~ X\'/'' 



is actually valid under the weaker assumptions of the proposition. It is then 
clear that Pa minimizes (|22p and that the other minimizers must satisfy \H — 
A 1 1/2 ('y — p^) — 0, which is the same as saying that the range of 7 — Pa is 
included in the kernel of iJ — A. D 

2.3.5. A representability criterion 

The aim of representability criteria is to identify sets of densities p that arise 
from admissible density matrices. For finite systems, if 76 ©iHiS, 0<7< 1, 
and Tr (— A7) < cxj, then p-^ > and ^/p^ <E H^ (R'') by the Hoffmann-Ostenhof 
inequality. Lieb's representability theorem [20|, Theorem 1.2] shows that these 
conditions are sufficient for a function p to be representable. 

In the ergodic case, we know that a density p must satisfy p > 0, ,/p G 

Hi and p <E Lg , by the Lieb-Thirring inequality ([7]). Clearly a stationary 
function p such that p > and ^/p e H^ is not necessarily the density of an 
ergodic density matrix with finite kinetic energy, since in general 

{p>0\VpeHl}<^L^. 

It is an interesting open problem to determine the exact representability condi- 
tions in the ergodic case. Theorem 12.131 below gives sufficient conditions for p 
to be representable. These conditions are also necessary for d = 1. 

Theorem 2.13 (A sufficient condition for representability). We assume that 
d > I. Let p be a function satisfying 

p > 0, p e LI and y/p e HI . 

Then, there exists a self-adjoint operator 7 in &^ 1 H 5, satisfying < 7 < 1 
and p^ = p a.s. 

Theorem 12.131 is proved in the Appendix, following ideas of Lieb [20| . 
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3. Yukawa and Coulomb interaction 

This section is devoted to the definition of the potential energy per unit 
volume of a stationary charge distribution /. In our setting, / will be p^ — /i, 
where fi is the nuclear charge distribution and pj the density associated with 
an electronic state 7. We will consider two types of interactions, namely the 
(long-range) Coulomb and the (short-range) Yukawa interactions. 

In dimension d > 1, the Coulomb self-interaction of a charge density / is 
given by 

D{fJ)^ I V{x)f{x)dx. 



where V is the Coulomb potential induced by / itself, which is solution to 
Poisson's equation 



/w ^ \s''-'\ f. 



(23) 



Here |5'''~^| is the Lebesgue measure of the unit sphere 5'''"^ (l'5'"| = 2, |5^| — 
27r, \S'^\ = Att). For later purposes, it is convenient to regularize this equation 
by adding a small mass m as follows : 



{-A + m')V=\S''-'\,f. 



(24) 



Whenever ttt, = or 771 > 0, we have the following formulas for the Coulomb 
(to = 0) and Yukawa (m > 0) self-energies: 



D^,{fJ) = \S' 



d-ll 



fiK) 



\K[ 



■dk 



^d-ll 



(-A- 



L2 



(25) 



%nix-y)f{x)fiy)dxdy (26) 



W„i{x-y)f{y)dy 



dx. (27) 



Here / is the Fourier transfornij of /. Of course we need appropriate decay and 
integrability assumptions on / to make the previous formulas meaningful. The 
Yukawa and Coulomb kernels are given by 



- m I x I 



Y^,ix) 



Ko{m\x\) and Yo{x) 

\^\~lp~m\x\ 



-\x\ ifd==l, 

-log(|x|) ifd = 2, 
|x|~i lid = 3, 



withA'o(r)=/o" 



sht 



dt the modified Bessel function of the second type [22 



The Coulomb potential is nothing but the limit of the Yukawa potential when 
the parameter to goes to 0. Similarly, the function Wm is defined by its Fourier 
transform 



WraiK) 






^In the whole paper we use the convention / {K) = {2tt) 2 J^^ f (x) 



"dx. 
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Using the integral representation x ^^'^ = 27r -'^ /o°°(-'^ + '*^) ^'^*' ^^ ^^^ ^'^^^ 



2 



W„,{x) = ^j=== I Y^^r+^{x)ds. (28) 



This can be used to compute Wm in some cases, or to simply deduce that, when 
771 > 0, Wm is positive, decays exponentially at infinity, and behaves at zero like 
\x\~'^ in dimension rf = 3, like |a:|~^ when d = 2 and like log |x| for d= 1. 

Our goal in this section is to define the Yukawa and Coulomb energies per 
unit volume for a stationary charge distribution /. Formally, this is just 



where V solves ([?5|) for m = or (|M|) for m > 0. We are implicitly using here the 
fact that the potential V is stationary when / has this property. Unfortunately, 
giving a meaning to Poisson's equation (P^ in the stochastic setting in not an 
easy task. Already when / is periodic, we know that this equation can only 
have a solution when /„ / = 0. Here the situation is even worse, as we explain 
below. To simplify matters, we first introduce the Yukawa energy per unit 
volume Dm (/, /) for to > and then we define the Coulomb energy as the limit 
of Dm{f, /) as 771 — > 0, when it exists. Thus we start by giving a clear meaning 
to the three possible formulas ([25]). (pS)) and (P7)) in the Yukawa case tti > 0. 
In the next section we introduce the stationary Laplacian — A^ which allows to 
write a formula similar to (1251) . 



S.l. The stationary Laplacian 

In this section we define an operator which we call the stationary Laplacian, 
which is nothing but the usual Laplacian in the x variable acting on L^{fl x Q), 
with stationary boundary conditions at the boundary of Q. Surprisingly, this 
operator does not seem to have been considered before. 

Let Aq be the operator on L^ defined by 

D{Ao) = LlnL'^{n,C^{R'^)), 

Ao/(a;,x) = —Af{uj,x) a.s. and a.e., V/ e D{Ao), 

where A = J2i=i ^x refers to the usual Laplace operator on C'^{M.'^) w.r.t. the 
X variable. Using stationarity, we obtain 



(/,v4o.9)L^=E(y Vf{;x)-Vgi;x)dx], yf,geD{Ao). 

Thus, Aq is a symmetric, non-negative operator on L^ with dense domain D{Aq). 
We denote by — A^ its Friedrichs extension [35|, theorem 4.1.5, p. 115], and call 
the operator As the stationary Laplacian. The form domain of the operator 
—As is Hi and its domain is H^. 

When il is finite, the spectrum of —As is purely discrete. If the probability 
space is defined as in Example 12.11 then tT(— As) = [0,-|-oo). Thanks to the 
ergodicity of the group action, one can prove that ker(— As) = span{l}. In 
contrast to the periodic case, there is (in general) no gap in the spectrum of 
—As above 0. In other words, there is no Poincare-Wirtinger type inequality 



18 



in Hi . This can be seen, for instance, by considering the sequence of functions 
^n{^,x) = n~'^/2$ (w,n~^a;), where $ {uj,x) = Y.kez<i ^('^fc(w))x(a; - k) with 
Y G L^ (il) and x G C^ (R'') with support in the unit cube Q and such that 
lo x(^) '^^ ^ '^^ These functions are such that ||$,i||^2 = 1 and ]E(/q $ti) = 
for any n g N, and j| V$„||(-^2id — > as n — > oo. 

That there is no Poincare-Wirtinger inequahty means that solving Poisson's 
equation ([^5]) in the stochastic (ergodic) setting is comphcated. Contrarily 
to the periodic case, it is not sufficient to ask that / <^ ker(— As), that is 
£(/„ /) = 0. If we are given f £ L^, then we see that there exists V £ L^ such 

that —AsV = \S'^^^\f if and only if / belongs to the range of — A^. In the next 
section we consider the simpler Yukawa equation (j24p . 

3.2. The Yukawa interaction 

Let m > 0. If / e Lg, we can define by analogy with ([^S]) 

D,MJ) = \S''-'\ {-A.,+mY\f\- (29) 



The operator (— A., + m^) ^ being bounded, D,„ is well-defined on L^. To set 
up our mean-field model for disordered crystals, we however need to extend the 
quadratic form Dm to a larger class of functions. Formal manipulations show 
that for a stationary function / 



Drn{fJ)=^{l I Y, 



Xx-y)f{-,x)f{-,y)dxdy 




Wm{x-y)f{-,y)dy 



dx . (30) 



The second formula is more suitable for a proper definition of -D,„. We claim 
that the function (W,„ * /) (w, x) := J^^ f (oj, y) Wm {x — y) dy is well-defined 
for all f € L\, and is in L\. This follows from the following elementary result. 



Lemma 3.1 (Convolution of stationary functions). Let J e L\{L'') and W € 
i[^^(]R'') such that 

for some 1 < p,q,t < oo. Then the function 

(W * f) (w, x):= [ f (w, y)W{x~ y) dy (31) 

belongs to Ll{L^) with 1 -|- 1/r = 1/p + 1/q, and 

WW * /||^,(^.) < C ||/||^.(^,) I ^ ||I^IL.(Q+,) 1 (32) 

\feez<i / 

for a constant C depending only on the dimension d. If 1 < p,q,r < oo, we can 
replace ||M^|lLP{n+fc) by the weak norm \\WtQ+k\\i^p in dSl])- 

19 



Proof. In order to prove the convergence of the integral in pip , we write 

/ \f{oo,y)\\W{x-y)\dy=Y. I \f {^,y)\\W {x - y)\dy 

= E / \f{rk{u:),y)\\W{x-y-k)\dy 

where we have used the stationarity of /. By the standard Young inequahty we 
have for a.e. x € Q and a.s. 



\fiTk{Lo),y)\\Wix^y~k)\dy 



Q 



LHQ) 



<l|W^(--fc)llL.(2Q)ll/(^fcM, 



li^Q) 



and therefore 

\f{Tk{Lo),y)\\W{x-y-k)\dy 
Q 



LiiL-) 



<l|W^(--fc)llLP(2Q)ll/llLi(L.). 



The rest foUows. The estimate with the weak norm ||M^lQ+fe||^p follows from 
the generalized Young inequality |29|. D 

Since Wm G L^(M.'^) when m > 0, Lemma [3TT] shows that Wm * f ^ L\ when 
f & L\. Now we can define 



D^ifJ):-- 



\W^ * /|- 



for any / in the space 



■DY:={fGLl\Wm*feLl} 



(33) 



(34) 



which we call the space of locally integrable functions with locally finite Yukawa 
energy. It is easy to see that the space Vy in fact does not depend on m > 0. 
It is a subspace of Lj , with associated norm ||/||^i +Di(/, fY^^, and a Banach 
space for this norm. 

Using Lemma 13.11 ((28)) and the known properties of W,„ , we deduce the 
following result. 

Corollary 3.2 (Some functions of Py). We have, in dimension d, 

(Ll{L^) ifd^l, 

LIdVyD IlI{L1), Vg>l ifd^2, 

[ 2.2(2.6/5) ,Jrf^3. 

When / S Lg we have that both Ym * f and Wm * f belong to L^, since 
Yjn and Wm are in L^{R'^) when m > 0. Thus we always have / € Vy and it 
is then an exercise to show that all the formulas for -D,„(/, /) in ^^ and ([50]) 
make sense and coincide. Indeed, we have 



Ym^f^\S'-'\{-A 



m 



f and Wm * I 



\S'^- 



2\-l/2 



-A, + ™2)-^'V. 
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3.3. The Coulomb interaction 

As mentioned previously, the Coulomb potential can be seen as the limit of 
the Yukawa potential when the parameter tti > goes to zero. More precisely, 
as < (-A,, + 7«2)-i/2 < (-A3 + 7n'2)"i/2 for all rn > in' > 0, the function 
m I— > Dm (/, /) is non-increasing on (0, +00), for any / £ 'Dy ■ It would therefore 
be natural to define the average Coulomb energy per unit volume as the limit 
of Dm (/, /) when m — )• 0, but we will proceed slightly differently. 

To simplify some later arguments, we define the Coulomb energy per unit 
volume by compensating the charge by a jellium background. This means we 
introduce for a stationary charge distribution / G Py 





D„{f,f):^ \im Dmlf- 

m->0 y 


-E 


X^) 


, ./ - E ( 


7^) 

Jq J 


J ' 




together w 


th the associated space 
















Vc = \l (^Vy\ lim An(/- 


-E 

v 


JQ J 


,./-IE 

V 


Ul 


)<oo}. 



of the locally integrable stationary charge distributions / with locally finite 
Coulomb energy (when compensated by a jellium background). We again em- 
phasize that Dc C Ll by construction. 

When f G L1 O Vy, the limit is finite if and only if / — ^{Jn f) belongs to 
the quadratic form domain of (— A^)^^, and we have by the functional calculus 



DoifJ)^ lim A, 

m— >-0 



S" 



d-ll 



./-: 



(-A.) 



/ , / - E / ./ 



/-E( / / 
Q 



Li 



For / only in Vc, the family (-A, + m^)-^/^{f - E{Jq f)) is Cauchy in L^ 
when m goes to zero and we still denote its limit by (— As)^^/^(/ — E(J„ /)). 

The following result means, in particular, that in the physically relevant case 
d = 3, a stationary function / € Ll{L^/^) whose charge and dipole moment in 
the unit cell Q vanish a.s., has a finite average Coulomb energy per unit volume. 

Proposition 3.3 (Some functions in I?c)- Let d < i and f he a function of 
Lg {L'i), with 9=1 if d = 1, q > 1 if d ~ 2 and 9=1 if d ^ 3, such that 



q{uj) 



f{uj,x)dx = and p{ui) 



xf{uj,x)dx = a.s. (35) 



Then, f € T>c. 



Proof. For the sake of brevity, we only detail the proof for d = 3. Let / be a 
function of L^ (i^) satisfying ([55)1 . As E(J„ /) = 0, we have for all to > 0, 



Dmlf-^i / / , / 



/ 



\x-y\ 



fi-,y)fi-,x)dydx 



<C\\f\\ 



LliL^) 



■E 
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where 

!• ^~7n\k+y-x\ 

An.fe(w) = / 7^— rf {uj,y) f (Tkiuj) ,x) dydx. 

jQxQ {fi + y - X\ 

Noticing that for all m > 0, {x,y) G Q x Q, and k £ Z'^ such that |A:| > 3, 



e 



-m\k+v-x\ __ g-m|fc| ^ ^e-'n|fc|(|fc + y - x] - |fc|) < m2e-™l'=l/^ 

and using the fact that q{uj) ~ a.s. we obtain that 

A„,,k{oj) = (1 + m|fc|)e^"l'=lB„,,fe(a;) + a„,fe(w) a.s., 
with 

""-^'^^^^L, \k + y~.\ '''' 

and 

|C™,.H|<2m2|fc|-ie-"l'=l/2||/(..,.)|UHQ)ll/(^^H'-)llLi(Q). 
It then follows from (1551) that 



Sm.fc(w) = / f{uj,x)f{Tk{uj),y)F{k,x~y)dxdy a.s., 

jQxQ 

where F(fc, /i) = ^ - (^ - £|^ + MSi^l!^) and e^ == k/\k\. Thanks to 
the multipole expansion formula (see e.g. |l9l . Lemma 9]), there exists a constant 
C such that for ah fc e M^ \ {0} and /i € K^ with fc + /i 7^0, \F{k,h)\ < TjMrr^:^- 
Therefore, 

|S™,;.MI < C\k\-' 11/ (^, OILifQ) 11/ (r, H , .)JL.(Q) a.s. 
Consequently, 



\fceZ3,|fc|>3 / \A:GZ3\{0} ' ' kel,^\{ 



-7n\k\/2 



M 
\{0} 

from which we infer that 

^g-m|fc|/2 



z.„(/-e(//),/-e(/^;))<c||/||^^,^,(i + J: 



2' 

m - 



k\ 

feGZ3\{0} ' ' 



for a constant C independent of /. As 



E-m|fc|/2 r g-\x\/2 

m — = / — -- — dx < 00, 
"'^u' \k\ /b3 \x\ 

we finally obtain that / e Pc- D 

The proof of Proposition 13.31 can be adapted to show that Do{f,f) is the 
limit of the supercell Coulomb energy per unit volume (see Section [5] and |6|), 
for any fixed / satisfying the neutrality assumptions psp . It is an open problem 
to prove the same for the functions / G Vc which do not satisfy psp . 
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3.4- Dual characterization 

The purpose of this section is to provide a useful characterization of the 
Yukawa and Coulomb spaces Vy and Vc by duality. Let us introduce the 
spaces of test functions 



and 



Ey = span {$x,y with F e L°°(f7), x e ^ (M'') } 



Ec = span{$x,y with Y eL°°{fl), x ^ So (K'^)} 



where S (R'^) is the Schwartz space, Sq (M'^) = {x e 5(M'^) | x e C;?°(E'* \ {0})}, 
and 

(f>Tj..y (w, a;) = > Y{Tk{uj))x{x — k), a.s. and a.e. 

The following says that Ey (resp. Ec) are dense in L^ (resp. in LP n 
ker(-A,)^). 

Lemma 3.4 (Density of _Ey and Ec)- For any 1 < p < oo, the set Ey is dense 
in LP and the set Ec is dense in LP ^ If e LP \E [/„ / j = [. 

Proof. We sketch here the proof of the density of Ec in if, and refer the reader 
to [III for further details. Let ip e (Lf)' = Zf', where p' = (1 - p^^)-^ is the 
conjugate exponent of p, be such that 



vr e L°°(f}), Vx e 5o(M^), e(I^<^^,y] ^0. 



(36) 



For Y e £°°(1^), we denote fyix) = E(Y{-)ip{-,x)). The function /y is in 
L^jjjf(R''), hence it is a tempered distribution: /y e 5'(M'^). In view of ([36| . 
we have for ah x S ^o, (^~^(/y), x)5'(R''),5(R'i) = 0. Therefore J^'^ify) is 
supported in {0}, which implies that 

\a\<N 

with A^ e N and c^ e C. It follows that 

fvix) = E Caa;", 

|q|<7V 

with Cq G C. As /y is in iunif(R''), all the coefficients Ca are equal to zero, 
except possibly cq, and /y is a constant. As y i— > /y is a continuous linear form 
on LP{n), there exists Z e LP (fl) such that for all Y e L°°(f7): E (yZ) = /y. 
It follows that for all x G K'', (p{uj,x) = -^(w). We know that any stationary 
function independent of a; is a.s. and a.e. constant |28| . As K{J„(f{uj,x)) ~ 0, 

we conclude that (p ~ 0, which proves that Ec is dense in L^. D 

It can be verified [18| that 

V$x,yG^y, i-^s + l)'-''5?x,Y^'^(-A+i)-^/-x,Y (37) 
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and, similarly, that 

V$x,y e Ec, (-A,)"5 ^^^ ^ $(-a)-i/^x,i'- (38) 

A straightforward consequence of Lemma 15^ and ([57)) -(|55 | is the following 
Corollary 3.5 (Dual characterization of 2?y and Vc)- Let f ^ L\. 

(i) // (—As + 1)^^"/, seen as a linear form on Ey , is continuous on 

{EvAl-lDLi, thenfeVy and Dra{fJ) = \S''-'\\\hA, + m^y^/^ f\\l^. 

(ii) IfK{jQf) = and {—As)^^''^f, seen as a linear form on Ec, is continu- 

ous on (EcW-Wli), then f e Vc and Doif, f) ^ \S''~A\\{-As)-'/y\\l,^. 

4. Stationary reduced Hartree-Fock model 

After these long preliminaries, we now introduce and study a reduced Hartree- 
Fock (rHF) model for crystals with nuclear charges randomly distributed fol- 
lowing a stationary function /i > 0. We typically think of /z being of the form 

/^(w, a;) = ^ qk{uj)xix - k- rjuiuj)) 



with J X = 1 and which describes a lattice of nuclei whose charges and positions 
are perturbed in an i.i.d. ergodic fashion. However in this work we do not 
want to restrict ourselves to /i's of this very specific form and for us /i is any 
non- negative stationary function in L\ . Our only restriction in this work is that 
we do not allow point-like charges. 

In Section 14.11 we define the minimization sets and the rHF energy func- 
tionals associated with the Yukawa interaction of parameter tti > on the one 
hand, and with the Coulomb interaction on the other hand. In Section 14.21 we 
prove the existence of a ground state density matrix 7, and the uniqueness of 
the associated ground state density p^. We then show in Section H751 that the 
m- Yukawa rHF model converges to the Coulomb rHF model when the param- 
eter 771 goes to 0. Finally, we prove in Section W^ that, in the Yukawa setting, 
any rHF ground state satisfies a self-consistent equation. 

In Section [SJ we will prove that, still in the Yukawa setting, the rHF model 
for disordered crystals we have introduced is in fact the thermodynamic limit 
of the supercell model. 

4..I. Presentation of the model 

As in the usual rHF model for perfect crystals (6|, the rHF model we propose 
consists in minimizing, on the set of admissible density matrices, an energy 
functional composed of two terms: the kinetic energy per unit volume and 
the average Coulomb (or Yukawa) energy per unit volume. This leads us to 
introduce the family of energy functionals 



^^,"1(7) =2— (^^'^)+ 2^"^^'^^'"''''^"'^^ 
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(39) 



with m = for Coulomb and to > for Yukawa. The sets of admissible density 
matrices are defined by 

IC^^Y = lie&,^,nS, 0<7<la.s., TLh)=^(f A , P^ - ^i e Vyl 

(40) 
in the Yukawa setting, and by 

/C^,C = {7 e ©1,1 n 5, 0<7<la.s., Tr(7)=EU pA , p^ - fi e vA 

in the Coulomb setting. The constraint Ir (7) = E(/q p) (neutrality condition) 
must be added in the latter setting since the average Coulomb energy per unit 
volume of a non globally neutral stationary charge distribution is infinite (recall 
that in our definition of Dq, we have added a jellium background to enforce 
the neutrality condition) . We also impose this constraint in the Yukawa setting 
for consistency. In our model it is not essential that ^ > but we keep this 
constraint for obvious physical reasons. 

The following lemma gives sufficient conditions on /^ > for the sets /C^i.f 
and /Cp.c to be non empty. 

Lemma 4.1. // /i G 2?y, then K,^y is non empty. // /i > satisfies the 
following conditions 

(i) p e Ll{L^), 

(ii) there exists e > such that \p{uj)\ < q{uj)(^^ — e) a.s., where q{uj) — 
/„ p{oj, x) dx and p (lu) = J„ xp{uj, x) dx, 

then /C^,c is non empty. 

Loosely speaking, the interpretation of the condition \p (w)| < q (w) (5 — e) 
is that the nuclei do not touch the boundary of Q too often. 

Proof. Let p G Vy and p := E(/q p) a.s. and a.e. It is clear that there exists 
a self-adjoint operator 7 G &_^ ^ such that < 7 < 1 a.s. and p^ — p. We can 
take for instance a free electron gas with constant density p, that is. 



7 = l(_,e](-A), with £=(fg^) p''' 



\s<^-~ 

This state is obviously ergodic since it is fully translation- invariant. Moreover 
it satisfies 

Tr (-A7) = ^ r^^i^ V'%-V. 

Besides, p ~ p £ Vy and therefore 7 G JC^,y. 

Suppose now that p satisfies conditions (i) and (ii) of the statement. Let p 
be the stationary function defined on Q by 

'0 ifg(a;) = 

p{oj,x)={ ,. f^-sl^\'^ (41) 

^&x{-ri^] otherwise. 
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Here d{uj) = dist{p (lu) / q (uj) , dQ) and x is any non-negative radial function 
of C^(R'^) with support in B(0,l/2), such that J^x^ = 1- We check that 
p ^ L^ {L^) n Lg. where q satisfy the conditions in Proposition 13.31 and ^ £ 
Hi. Therefore, by the representabihty Theorem l2.131 there exists a self-adjoint 
operator 7 £ &_^ ^ such that < 7 < 1 a.s. and p^ = p. Moreover, 
J„{p{iu,x) — p{u!,x)) dx = and J„x{p{uj,x) — p{uj,x)) dx = 0. It follows 
from Proposition 13.31 that p — fi E Vc, and therefore that 7 e ^ti,c- D 

4.. 2. Existence of a ground state 

Now that we have properly defined the rHF energy, it is natural to look for 
ground states, that is, niinimizers of £^^m on IC^y/c- The ground state energy 
of a disordered crystal is defined by 



Ifi.ni = inf {£^,™(7), 7 e JC^x) (42) 



with 771 > 0, in the Yukawa case, and by 



J^,o = inf{fp,o(7), TG-'C^^c} (43) 



in the Coulomb case. 

Theorem 4.2 (Existence of ergodic ground states). Let < p E Ll. If ICf^,Y 
(resp. K,^^c) is non empty, then |7^ (resp. \4-^ ) has a minimizer and all the 
minimizers share the same density. 

The proof of Theorem 14.21 is based on the weak-compactness of ICfj^^/c 
(Proposition 12. 9|) . and on the characterization of the spaces Vc/y by dual- 
ity f Corollary 13. 5p . We recall that in Lemma BTT] above, we have given natural 
conditions which guarantee that /C^.y/c is non empty. 

Proof. Let m > Q and let (7^) be a minimizing sequence for /^^m- As the 
functional f^^m is the sum of two non-negative terms, these two terms must 
be uniformly bounded. Since Tr(— A7„) and 2x(7n) — ^ijn t^) ^^^ bounded, 
we can apply Proposition 12.91 and extract a subsequence (denoted the same 
for simplicity), such that 7„ ^* 7, with all the convergence properties of the 
statement of Proposition [^7^ In particular, we have 

Tt (-A7) < lim infix (-^Jn) and II (7) = E ( / /i 

Similarly, we know that z„ :~ Wm * {p-y„ ^ A*) = (^As + m'^)^^i'^{p^^ — /i) is a 
bounded sequence in L^. Thus we can extract another subsequence such that 
z„ ^ 2; weakly in L^. 

Passing to weak limits using that p^^ -^ p^ in Ls , it is readily checked 
that for any $ € E^jy 

p^-/^, (-A, +m2)"^$\ ^ =(z,$)^,. 

Therefore, using the lower semi-continuity of the L^-norm, we obtain 

m'^y~^{p^-p) = ||z||^2 <liminf (-A, + m2)"^(p-y„ - m) 
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We deduce from Corollary 13.51 that p^ — i^i E 'Dq/y a-nd that 
■Sp.mlT) < liminf £^,.,„(7„) = /^,„. 

n— ^oo 

Thus, 7 is a minimizer of (|^^ (resp. (|^5)) ). 

Let us now prove the uniqueness of the minimizing density p^ . Assume that 
7i and 72 are two minimizers of (j42p (resp. (j43p ). A simple calculation shows 
that 

= J/j,!?! ~ "^^m [Pji ~ P72 1 P71 ~ P72J ■ 

As /^,„i is the infimum of £fj,^„i and as (71 + 72)/2 belongs to the minimization 
set JC^^c/Yi we deduce that ||(— As + m^)^2 (p^^ — p^^) ||^2 = 0. Thus (— A^ + 
rri^T^ (P71 - P72) = 0- For all $ G i^c, (-A^ + m2)i/2$ g £:^ and 

((-A, + m'y'' {p^, - p^J , (-A, + w?)^-<^)^^_ = 0. 

Hence, E(/q(p^i - ^72)$) = for aU $ e E:^- As i^c is dense in l\'^'^''^ n {1}-^ 
(see Lemma 133]) and as, in addition, E(/q(p^j ~ Pn)) — Oj ^^ conclude that 

P71 = P72- n 

/^.5. From Yukawa to Coulomb 

In this section, we prove that the ground state energy of the Yukawa problem 
converges to the ground state energy of the Coulomb problem as the parameter 
m goes to 0. The result essentially follows from our definition of the Coulomb 
energy Dq as the limit of Dm when m ^- 0. 

Theorem 4.3 (Convergence of Yukawa to Coulomb). Let < p E Ll be such 
that /C^.c 7^ 0. The function m i— > /^,m is decreasing and continuous on [0, +00). 
In particular, we have 

lim Iu,m = In.o- 



Moreover, if for each m > 0, jrn is a minimizer of {4'/^ , then the family (7m)„i>o 
converges, up to extraction, to some minimizer 70 of C^, in the same fashion 
as in Provosition \2.9[ 



Proof. That m 1— >■ /^^^ is decreasing and continuous on (0, +00) is easy to check 
(the strict monotonicity follows from the existence of minimizers). For / € Dc 
such that E(/q /) = 0, we have Dm{f, f) < -Do(/, /) for all m > 0. It follows 
that 

V7 e /C,j,c, Vm > 0, fM,m(7) < '^^,0(7) < 00, , 

and therefore that 

If.,rn < /m,0- (44) 

This proves that limm_>.o+ I^^m < If^fi- 

For 771 > 0, we denote by 7^ a minimizer of (|^^ . We deduce from pi)) that 
there exists a positive constant C such that, for all m > 0, Tr (— A7„i) < C and 
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||(— As + m^) ^^'^{p-f,„ ~ m)|| r2 ^ C". Reasoning as in the proof of Theorem 14. 2 [ 
we can extract a subsequence {'ymk)keN with rrik \ 0, such that there exists 
J G K. with 



IL (7) - IImIIl;, S(-A7) < HminfS (-A7mJ , 

k:— >-oo 



and 



(-As)-Mp7-a^) 



< lim inf 

L^ k~>oo 



i^As + ml) 3 (^p^^^ -^) 



i? 



This proves that 7 G /C^j.c a-i^d that 



Iii.Q < £11.0(1) < hminffp.,„j7„J = lim /^,,„ < /^,o, 

/c— >-oo Tn— >-0 



which concludes the proof of the theorem. 



□ 



4-.4- Self- consistent field equation 

In this section, we define the mean-field Hamiltonian H = — ^A + V as- 
sociated with the ground state for tti > (Yukawa interaction), and we prove 
that any ground state of (|42|l satisfies a self-consistent field equation. The same 
holds formally in the Coulomb case but, unfortunately, we are not able to give 
a rigorous meaning to the Coulomb potential V. For this reason we consider a 
fixed parameter ?7i > in the rest of the section. 

We introduce the stationary mean-field potential V defined by 



V{uj,x) 



Ymix - y) (p^„ - ii) (w, y) dy, 



(45) 



where pj^^ is the common density of the minimizers of (j42p . The following says 
that, under the appropriate assumptions on ^, y is a well-defined stationary 
function such that the associated random Schrodinger operator H ~ — ^A -I- V 
is also well defined. 

Lemma 4.4 (Mean-field random Schrodinger operator). Let d G {1, 2, 3}, m > 

and < fj. E Ls D Dy ■ Let p^^ he the (unique) ground state electronic 

density for the Yukawa minimization problem (|42p . obtained in Theorem\ 
and V the associated mean- field potential defined in (|45p . Then we have 



V G { 



Ll{L°°) ford=l, 

Ll{L°-) ford = 2, 

_L^^^(L°°)nL2(L6) ford =^3, 



(46) 



and the random Schrodinger operator H := —^A + V is almost surely essentially 
self-adjoint on C^{M.'^). In dimension d = 3, if p, £ LJ {L^) j then we also 



V^ e Ll/'^. 



(47) 



Let us emphasize that H is a, uniquely defined operator since p^^ is itself 



unique. Note that under the sole assumption that n € L 



1+2/d 



in dimensions 



d = 1,2 we have p. e Vy by CoroUarv 13.21 In dimension d = 3, the additional 
hypothesis p G Lg' (L^) ensures that p <E Vy, by CoroUarv 13.21 and the fact 
thatLf/'(Li)nLs'/'cL2(L6/5). 
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Proof. As we know that p^^^ £ Lg , P5)l and (|47p follow from Lemma [3.11 
and the fact that V = cWm * {Wm * (P7„ — m)) with IVm * [p-tm. ^ m) "= ^s 
since p^,^ — fJ^ & 2?y . We know from [7|, Proposition V.3.2, p. 258] that —5 A + V 
is essentially self-adjoint on C^(R'^) when y g LI{LP) for some p > 2 and 
r > dp/{2{p — 2)). In our case we can apply this with (p, r) = (3, 3) for rf = 1, 
(p, r) = (5, 2) for d = 2 and (p, r) = (21, 5/3) for d = 3. D 

The following now gives the self-consistent equation satisfied by a minimizer 

7m- 

Proposition 4.5 (Self-consistent equation). Let d G {1,2,3}, to > and < 

fj, & Ls . Suppose also that p, € Lg (L^) if d = 3. T/ien there exists ep G M, 
called the Fermi level, such that any minimizer 7,„ of the Yukawa minimization 
problem (j42p is of the form 

7m = 1(-oo,£f)(^) +'5' 

/or some ergodic self-adjoint operator S satisfying < 5 < Ifep} (^)- 

Since H is uniquely defined, we deduce that two different minimizers need to 
have different operators S's at the Fermi level ep. In particular, when e-p is not 
an eigenvalue of H, we deduce that "f„i = l(-oo,eF)(-^) i^ ^^^^ unique minimizer 
of dH. 

Proof. As fi (z 'Dy, (|42p has a minimizer 7 by Theorem l4.2l The Euler inequality 
associated with the convex optimization problem P^ then reads: 

Vy e /C^,y, is (-A(7' - 7)) + D„.{py - p^, p^-p)> 0. 

For g e R+, we set 

E{q)= inf flTr{-A{y-^)) + D„,{py-p^,p^-fi)]. 

Tr(7')=<? 

It is easily checked that the function E is convex on R_|-, hence left and right 
differentiable everywhere. Also, for any 



ep e 



E' [Ei / n] -0] ,E' {e{ / p] +0 



(48) 



where E'{E{J„p) — 0) and E'{]&{Jqp) + 0) respectively denote the left limit 
and the right limit of the non-decreasing function E' at E(/q p), we have 

ITL (-A(7' - 7)) + Drnipy - P7. P7 - A^) - epll (Y - 7) > 

for any ergodic operator 7' e /C such that p^r E Vy ■ As p-^ E Vy , VJj = y,„ * p G 
is , and Pj' e is for any 7' e /C, the above inequality actually holds 

for any 7' G /C. In addition. 



-Cm(P7' -P7,P7 -M) =E / y(p7'-p7) 
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in R+ U {+cxj}. Taking now 7' = l(_oo,eF) (-^)' which belongs to /C by Proposi- 
tion [nHH and using Proposition 12. 121 leads to 

< ^II (-A(7' - 7)) + An(Pr - Pi^Pi - m) - £F^(7' - 7) 



Hence, 7 = 7' + (5 with 5 as in the statement. D 

The following result is a consequence of Proposition 14 . 5 1 and of the Feynman- 
Kac formula. 

Corollary 4.6. If fj, £ L°°{p.xR'^), then, for each m > 0, the common 
density p„i of the niinimizers of the Yukawa minimization problem (j42p is in 

5. Thermodynamic limit in the Yukawa case 

The purpose of this section is to provide a mathematical justification of the 
Yukawa model (|42p by means of a thermodynamic limit. So far, we did not 
manage to extend the results below to the Coulomb case. 

Let us quickly recall that the thermodynamic limit problem consists in study- 
ing the behavior of the energy per unit volume (as well as, possibly, the ground 
state itself and some other properties like the mean-field potential, etc) when 
the system is confined to a box with chosen boundary conditions and when the 
size of the box is increased towards infinity. 

For a perfect (unperturbed) crystal, the existence of the limit in the many- 
body case goes back to Fefferman [13J, after the fundamental work of Lieb and 
Lebowitz [2l|. A new proof of this recently appeared in |15| . However, for 
the many-body Schrodinger equation, the value of the limiting energy per unit 
volume is unknown. For effective theories like of Thomas-Fermi or Hartree- 
Fock type, it is often possible to identify the limit and to prove the convergence 



of ground states. In [25|, Lieb and Simon prove that, for the Thomas-Fermi 
model, the energy per unit volume and the ground state density of a perfect 
crystal are obtained by solving a certain periodic Thomas-Fermi model on the 
unit cell of the crystal. The same conclusion has been reached by Catto, Le Bris 
and Lions for the Thomas- Fermi- von Weizsacker model |8| , and for the reduced 
Hartree-Fock (rHF) model [9| we focus on in the present work. 

In the stochastic case, Veniaminov has initiated in J37| the study of the 
thermodynamic limit of random quantum systems, but with short range inter- 
actions. The case of a random Coulomb crystal was recently tackled by Blanc 
and the third author of this article in [3|. Blanc, Le Bris and Lions had already 
considered Thomas- Fermi like models in [2|, for which they could also identify 
the limit. 

Here we follow [6| and we consider the so-called supercell model. We put the 
system in a box F^ — [— L/2, L/2) of side L G N\ {0}, with periodic boundary 
conditions. When to > 0, we show that the ground states converge, when L 
goes to infinity, to a ground state of problem (^^ (up to extraction and in a 
sense that will be made precise later). 
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Let 771 > be fixed for the rest of the section. We mtroduce the Hilbert 
space 

ipcr (ri) = {ve iL (K') I V (iZ)'* -periodic} . 
The Fourier coefficients of a function / e L^^^. (Tl) are defined by 

We denote by — A^ and Pj.L^ ^ 1^ i ^ d, the self-adjoint operators on Lii^^ (F^) 
defined by 

4 (-Ai/) = I Af 4 (/) , and 4 (P,, i/) = fc,4 (/) , VA' e (^^z' 

For fc S Z'^. we denote as before by Uk the translation operators on L^^^ (M'') 
defined by Ukf {x) — f {x + k). For any /, g E Lii^^ (^l)i we set 

^™.L(/,5) = |^''-^|((-Ai + m2)-V,(-AL+m')"^9)^^ ^^ ^ (49) 

= E ,if , 4(7y4(g)- / f Y,n{x~y)f{x)g{y)dxdy. 

Denoting by ©i.l (resp. Sl) the space of the trace class (resp. bounded self- 
adjoint) operators on L'i^^.{rL), the set of admissible electronic states for the 
supercell model is then 

ICL = {lLe6ixnSL, 0<7<1, Tri2^^(r^)(-AL7i)<c»}. 

For any w G fi, we denote by ^l{^, ') the (AZ) -periodic nuclear distribution 
which is equal to /i(w, •) on F^, and by (5^,„ the (w-dependent) energy functional 
defined on K,l by 

^il^mi^^L) = -Tri2_^^(r^) (-Ai,7L) + -£'m,L(p7i, -Mi(w,-), P71, -ML(t^,-))- 

Let £f be as in Proposition 14.51 For any cu E il, the ground state energy of the 
system in the box of size A with Fermi level sp is given by 






(cj) = inf |£^^,„(a;,7L) - epTris^^^jr^) {jl) , 1l e /Cl| 



(50) 



Proposition 5.1. Let jj. E L^. For each A G N \ {0}, ([5D|) /las a minimizer, 
and all the minimizers of (j50p share the same density. 



Proof. The proof follows the same lines as the proof of [9|, Theorem 2.1], re- 
placing the periodic Coulomb kernel by the periodic Yukawa kernel Ym,L (x) = 

2_/fce(LZ)'' ^™(^ ~ ^)- '-' 
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On the other hand, the ground state energy of the fuh space ergodic problem 
with Fermi level ep is defined by 

-^M.m.EF = inf {£^l,rni1) " ^fTl il) , 7 € K-y} , (51) 

where £fj_,m, is given by (j39p . 

ICy ■■= {7 e ©1,1 n 5, < 7 < 1 a.s., p-y- fie Vy] (52) 

(the neutrality constraint has been removed compared to /C^i.y defined before 
in pp]) ). It is a classical result of convex optimization that (|42p and ([5T|l have 
the same minimizers. 

Theorem 5.2 (Thermodynamic limit for m > 0). Let fi E L^. We have 









To prove Theorem 15.21 we first establish preliminary estimates in Proposi- 
tion [5T3l Then, we prove a lower bound in expectation in Proposition 15.41 and 
an almost sure upper bound in Proposition 15.61 We then conclude the proof of 
Theorem 15.21 using Lemma 15.71 

In order to adapt our proof to the Coulomb case, we would need some 
estimates on the Coulomb potential Vl in the box F^. It is reasonable to 
believe that screening effects will make {Vl) bounded in, say, L^{p.,Ll^^^^{W^)). 
For a very general arrangement of the nuclei, bounds of this type are known in 
Thomas-Fermi theory (see [2|, Theorem 7], which is taken from Brezis' paper [5|) 
and in Thomas- Fermi- von Weizsacker theory J8|, Theorem 6.10], but they have 
not yet been established in reduced Hartree-Fock theory. Proving such bounds 
is of considerable interest, but it is beyond the scope of this paper. 

Proposition 5.3 (Upper bounds). Let /i G L^ and let Jl{^) be a minimizer 
'^f -^u rn epi'^)- Then, there exists C > and a sequence of integrahle random 
variables (Zl) converging to some Z £ L^ (il) a.s. and in L^ (n), such that 

I^m^eM + D,n,L {Py, (w, ') ,^7. (^> ')) < C L" Z L (uj) a.S., (53) 

E{TTLl^^(^r,){l~^LhL)+^{Drn,L{Pyr,~f^L,Pyr,-pL)) < C L" (54) 

for allLen\{Q}. 

Proof. Taking 7^ = as a trial state, we obtain that, almost surely, 

^-^^^1^0^ < ^ A„.L {pL (C., •) , /^L (W, •)) < ^^L (^) , (55) 

where Z^ = L^'^ Jp p^ converges to ]E(/q p'^), a.s. and in L^ (51), by the ergodic 
theorem. Besides, we have for any a G K and any 7^ G /C^, 

Tri.^^(r,) ((-Ai - a)^^) > -Tri.^^.(r,)( - Al - a) _ > -CU' (56) 



where C may depend on a and d, but not on 7^. The bounds (j53p and i|54 
follow from (155p . (|55|) . the positivity of each term of £^„(a;, •), and the fact that 
E(Zl) = ||mIIl2 is independent of L. D 
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Proposition 5.4 (Lower bound in average). Let fj, £ Lj. Then 

hmmt —: > If^.m,eY- 

The following definition introduced in [2| will be used repeatedly in the proof 
of Proposition [?7 



Definition 5.5. For a function 17 : fi x M'^ — > C and L E N, we call the tilde- 
transform g of g the following function 

g{uj,x) = —^ y^ g {T-k {<-^) ,x + k) U.S. and a.e. (57) 

fcerinz<i 

We can now write the 
Proof of Proposition \5.4\ Let 7l(w) be a minimizer of (|50p and set 

fceri.nZ'' 

Notice that p^^ = p^^ where the latter is the tilde-transform defined in ([571) . 
For any L G N \ {0}, we define the operator 

where ip^ is the (L'Z) -periodic function equal to tp on F^. It is easily checked 
that 7^ is self-adjoint and that < 7^ < 1- Thus, the family (7^) is bounded 
in L°° {n,B). Up to extraction of a subsequence, there exists an operator 7 G 
L°° (ri,i3) such that 7^ converges weakly-* to 7. Moreover, 7 is self-adjoint and 
< 7 < 1 a.s. Besides, 7^ {u!,x,y) = 71, {uj,x,y) a.s. and a.e. on 51 x F/, x F^. 
In the following, we will show that 7 G /Cy and that £f^ „i (7) ~ epTr (7) < 
liminfL-'^E(/^,„,,^(.)). 

Step 1. The operator^ is ergodic. Arguing like in the proof of Proposition 12.91 
it is sufficient to show that for aU u e L^{Vt), ip,ij £ C^{R'^) and R G Z'', 

e(^ {^{TRiu)) - UR^{iu)U*j,)^,^j)^,) = 0. (58) 

Let u, If, ip and R as above and L E N. We have 

fce(ri,A(ri.-i-fl))nZ'i 
where AAB := {A\ B)U {B\ A). Hence, for L sufficiently large, we have 

\E{u{{y^{TR-)~URj'j^U*R)^,^)L2)\^\E{u{i^L{TR-)-UR^LU*R)^L,^L)LHTr,)] 

The left side converges to E (u((7 o tr~ Ur'jU'^) ip, iji) l'^): a-nd the right side 
decays as L~^ . Thus, ((55)) is proved. 
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Step 2. We have 

S(7) = lim ^ ^7^^ '-. (59) 

Thanks to the estimate ([M)) . for any x £ W^'°°{M.'^), there exists a constant C 
such that for all L € N \ {0}, we have 

d 

|E (Tr (xiLX))\ + E 1^ (Tr {xPj%PjX))\ < C, 

Following the proof of Proposition 12.91 we can show that 

E(«Tr ixix)) = lim E (u f pyx^) (60) 

for aU u e L°°{n) and all x e if (K'*)- Choosing it = 1 and x = Ig, we get 

S(7)= hm fif /p^^ 
Finally, we remark that 

which concludes the proof of ([5^ . 

Step 3. The sequence {p^j^) converges weakly to p^ in L^'^'^''^{il, L^^^ (R'')). 
By (|60p. we obtain 

L-i-oo \ Jgd / \^ Jgd 

for all u e L°°(J7) and all x <= C*^ (R''). To proceed as in the proof of Proposi- 
tionHH we only need to show that {p;y^ ) is bounded in L'^+'^Z'^ [n, L^+Vd (^Bi)) , 
independently of L, for any compact set Bj — Uks^iiQ + k), with / C Z'' such 
that card(/) < oo. This bound now follows from the convexity of the function 
X I— 7- x^'^'^''^ and from the Lieb-Thirring inequality in a box |14| 

Step 4- We have 

E(TrLiJ^ALjL) 
Tr (-A7) <liminf ^ 



L^oo L<^ 
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As 7^ converges weakly-* in L°° {^,B) to 7, we can argue like in the proof of 
Proposition 12.91 and get 

d 

S(-A7) < IminfEEz^ ^ ^{{U:^n,L,Pj,LlLPj,LU*,ipn,L)mQ)) 

j=i neN ker^nz'' 
1 '^ 
i=i 
= ^'i'iS^ Ld^ (TrL^^^(r^)(-Ai7L)) , 
where we have used that the operators Pj^l commute with the translations Uk- 
Step 5. We have 



E 



Dm{p-f -~ ^J■l P-f — aO ^ liminf ■ 






(62) 



We denote by /^ = p^^^ — fi^ and f ^ pj — p. It follows from a simple convexity 
argument that for all fc G Z'', 



E 



'-AL + m'] ^/, 



. < — 7IE 

L^(Q+k)l L'^ 



'-Ai + m^) Vi 



As 



(-Ai+m2) Vl =|^''"'| Dra.LifLjL), 



we obtain that for all fc e Z'', 



-AL+m^ V. 



L2(nx(Q+A;)) 



< *^j lE(An,L(/L,/L))<g. 



Therefore, there exists a function z £ L^(r2, L^jjjf(]R'')) such that, up to extrac- 
tion, (-Al + to^)"^/^/l converges weakly to 2; in L'^(Sl,Ll^^^{M.'^)). By the 
weak lower semi-continuity of the L^-norm, we have 



E I / z^\ < liminf E 



-Ai+m^) Vi 



LHQ) 






We are going to show that z ~ (— A^ + m?) ^1"^ (p^ — /x), which will conclude 
the proof. To do so. we just need to check that for any u G L^^'^l'^ (J7) and 
XGC-fM''), 



^lirn^E [u \ X (-Al + m^) ^ /^ | = E ( w 



-A + m^) 'x)/ 



(63) 

Let u and ^ be such functions. Reasoning as in Step 1, we notice that the 
tilde-transform /l^ converges weakly to /i in L^^'^I'^{^^L^^^ (R'*)). Then, we 
proceed in two steps. First, we show that 



X(-Ai+77Z2) Vl 



'?/l, 
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where rj = (— A+m^) ^^^x- Recall that, for any h G 5(M''), the function defined 
by hL (x) = EkeiLzr h{x-k) is in L^^, (r^) with c|^ (/i^) = {2iT/LY/^h [K). 
For L sufficiently large, we therefore have 



x(-Ai+m2) - h 



x)' Y. 



'{K)c'i,{fL) 



Ke^ 



Kr +m^ 



tY E ^n^y 



K 



Kei 



= J2 ''K(-nL)c'k 



K£i 






Next, using the fact that r\ g 5(R'^), the weak convergence of /^ to / in 

.1+2 



ii+2/''(^,Lj+^/'*(R'^)), and the bound (gl]), we obtain 



u / n!i 



u vf 



This concludes the proof of (|55)) . hence of (|5^ . 

Proposition 5.6 (Almost sure upper bound). Let /j. E L^. Then, 



hm sup -^ J 



(c) 



</, 



fi,rn,eF i ^•^- 



D 



(64) 



Proof. We will prove ((M|) assuming that /x S L°° (fi x R'') ; the generalization is 
obtained by an e/2 argument using (j53|) and ((56)) . Let 7 be a minimizer of ([ST]). 
By the ergodic theorem, there exists i7' C Vt, with P (17') ~ 1, such that on Vl' 



lim — 7 / p-., = E 
for any \ < j < d^ and 



P7 



L-i-oo L** 



PP,1P, 



E 



PP,7P, 



lim 



-A, 



,2\~2 



(P7 - /^) 



^"(ri.) 



Ld 



(-A, 



,2^-2 



(P7 ~ M) 



(65) 



Let Wo G ^' be fixed for the rest of the proof. Let < xi, < 1 be a sequence 
of localization functions of C^ (R'')^ which equals 1 on F^-i, has its support 
in Tl, and satisfies |Vxl| < C For L S N\ {0}, we introduce the operators 
jI : L"^ {W^) -^ I? (R'^) and 7^ : i^^, (Fl) -^ LI^^{Tl), whose kernels are 
given by 

'yLix,y) = XL{x)j{uJo,x,y)xL{y)andjL{x,y)= ^ jl{x + j,y + k). 

Using similar techniques to the ones used in the proof of Proposition 12.81 one 
can show that 

I^c.1/ L P-'^=^(fpX (66) 
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and that 

r^™ Td^rL^ (Fz.) (-Al7l) = Tt (-A7) . 

We now turn to the convergence of the potential energy, i.e. 
hm T-rAn.L {9l,9l) = Dm (/, /) , 



(67) 



(68) 



where f = p^ — fi and g^ = p-^^ — ^i. We introduce the auxiliary function f^, 
defined as the LZ'^-periodic function equal to / on r^. We first prove that 

r^™ Td i^^,L {9l,9l) - Dm.L UlJl)) = 0. 

Indeed, rewriting ql as ql = xi,por/i + (xi,per ^ IJA^i. with the definition 
XL.per = Efce(LZ)'* XL (' + fc), we have 



-Ai + m^) ^-{9l-Il) 



< 



L^Tr^) 



m ^||(5L-/L)lli2(r^) 



< ^ ^||(xi,pcr- 1) (/l-Ml)| 



i^=^(ri.) 



which is a o{L'^). Then, we prove that 

y™' Fd = -^™ y^^ /) ' 

To do so, in view of (|55|) it is sufficient to show that 



5* 



d-ll 



O-L 



Ld 



^AL + m'] - f, 



LHTr^) 



-A,+m^ -f 



LHTl) 



(69) 



(70) 



tends to zero. This follows from the fact that 
f 



tti = -^ / dxf (x) 



dyYm{x~y){fL{y)-I{y))., 



f e L°" {fix W^) and y„ e i^ (R'^). This completes the proof of dSH]). Com- 
bining ^^, (inZI) and (IHHl), we end up with 



tL ( \ 

hmsup -^ — — ^ < li 

L— i-oo ^ 



fj>,m,€F 



n 



for every wq G f^', which concludes the proof of Proposition 

We complete the proof of Theorem 15.21 using Lemma 15.71 below applied to 
Xl (w) = L-'^I^^.e^itu) and the bound (P)). 

Lemma 5.7. Let (X„)^gp^ &e a sequence of random variables in L^ (f2) and X G 
L^ (fl). We assume that there exists a sequence of random variables (.^n)jigN 
converging in L^ (fl) to Z E L^ {fl) such that 

• liminfE(X„) >E{X) 

n— >oo 

• limsupX„ < X a.s. 



37 



• ^n 1^ Z„ a.s. 

Then, Xn — > X strongly in L^ (O) as n —> oo. 

Proof. Replacing Xn by Xn — X , we can assume without loss of generality that 
X ^0. We then write X„ = (X„)+ - (Xn)-- We first notice that (X„)+ -> 
a.s. By the dominated convergence theorem with "moving bound" (see e.g. 
[2 2 . Theorem f .8]), we conclude that {Xn)+ — > in L^ (fi). By the liminf 
condition, we have limsup^^^^ E ((X„)_) < 0. As {Xn)- > 0, we conclude 
that (Xn)- -!> in Li (17). Finally, E(|X„|) = E((X„)+) + E ((X„)_) tends to 
0. D 



Appendix: Proof of Theorem 12.131 



Here we write the proof of Theorem l2.13l This transposition of Lieb's repre- 
sentability theorem to the ergodic setting claims that for any p satisfying p > 0, 
p £ Ll and y/p G H^, there exists a self-adjoint operator 7 G &_i i, such that 
< 7 < 1 and p^ ~ p. 



Proof of Theorem \2.13[ We start with the case d 
tions ipo,ipi £ C^{M.'^) satisfying 



1. We consider two func- 



• supp((po) C 



2' 21 



and supp((y9i) C [0, 1], 
• SfcGZ Vk = 2 where ip2k{-) = ^o{- - k) and (p2k+i{-) ^ (pi{- - k). 

We denote by 

Pk{LO,x) :^ p{LL>,x)ipk{x), 

and observe that p = X^feez P^l"^- Let Nk{uj) ~ L^ Pki^, x) dx. For each fc e Z, 
we set ipj^k = for all j e Z if Nk (uj) = 0, and 



ipj^k{,(^,x) 



\/pk{i^,x) 



exp 



2iTrj 



Pk {uj,t)dt 



otherwise. We then introduce the density matrix 

7fe = ^'rij^kW3,k){^j,k\^ 
j&i 

where nj,fc(a;) = lj<jVfc(cu) + {Nk{i^) - [iVfc(w)])lj=[Ar^(„)]+i. 

Each 7fe is in /C = {7 G ©^ ^ n 5 | < 7^ < 1 a.s.} and p.y^^{u}, •) = Pfc(w, •) 
a.s. As the supports of the kernels of jk and jk+21 are disjoints for all k,l G Z, 
the operators 7e ~ X)fcez72fc and 70 = X]feez72fe+i are in /C. By convexity, so 
is 7 = -^^Y^- ^^ is finally easily checked that pj = p. 

We now turn to the case d = 2. In the same spirit as for d = 1, we cover the 
space with a finite number of periodic patterns, in such a way that the elements 
of each pattern do not intersect (see Figure [T]). For example, let 



Ao 



r 5 5 \ 


2 


1 2\ 


1 1\ 


\ 1 1\ 


1 2\ 






, Bq = 


-,- X 
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— ,- X 


-,- 
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[ 12'12J 




3 3/ 


[ 4' Aj 


[ 4'4; 


3 3/ 





1 5 
6'6 
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The Z^-translations of these sets Ik ~ Iq + k, I E {A, B, C}, satisfy Ik Ij ~ ' 



for k ^ j and U^gz^^fc UBkUCk 



regular functions ((/sfjfcgz'^j I € {A, B, C}, such that 



Next, we consider three sequences of 



ifi > 0, suppifi) c Ik 



and J2^^ + <fi + ^^ ^3. 
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, iTT 


1 1 


1 ' ! 


1 n ~i 
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1 1 


1 1 


1 1 
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1 1 


1 












_i_ J 








::_Lr::::::: 



Figure 1: Covering in dimension d = 2 used in the proof of Theorem 12.131 

Repeating the argument detailed above in the one-dimensional case, we de- 
fine 7/, for / e {A, B, C}, and 7 ~ J2ieiA b c\ lil'^ and we check that p^ ~ p 
and that 7 satisfies the desired conditions. We proceed similarly for d > 3. D 
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